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We  observe  intense  pulses  of  far-infrared  electromagnetic  radiation  emitted  from  arrays  of  InAs  nanowires.  The 
terahertz  radiation  power  efficiency  of  these  structures  is  ~15  times  higher  than  a  planar  InAs  substrate.  This  is 
explained  by  the  preferential  orientation  of  coherent  plasma  motion  to  the  wire  surface,  which  overcomes  radiation 
trapping  by  total-internal  reflection.  We  present  evidence  that  this  radiation  originates  from  a  low-energy  acoustic 
surface  plasmon  mode  of  the  nanowire.  This  is  supported  by  independent  measurements  of  electronic  transport 
on  individual  nanowires,  ultrafast  terahertz  spectroscopy,  and  theoretical  analysis.  Our  combined  experiments 
and  analysis  further  indicate  that  these  plasmon  modes  are  specific  to  high  aspect  ratio  geometries. 
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I.  MOTIVATION 

Manipulation  of  electromagnetic  fields  at  subwavelength 
dimensions  is  an  emerging  trend  in  nanoscience  and  nan¬ 
otechnology.  This  is  realized  by  modification  and  control 
of  material  properties  and  geometries  on  the  nanoscale.1 
For  instance,  the  response  of  free  carriers  to  external  per¬ 
turbations  gives  rise  to  modified  behavior  of  slab  plasmon 
and  plasmon-polariton  modes  in  confined  geometries.2  The 
conversion  of  electromagnetic  waves  into  propagating  charge 
carrier  waves  allows  geometric  control,  localization,3  and 
guiding  of  light4  at  sub  wavelength  dimensions;  this  has 
provided  the  impetus  for  the  emerging  field  of  plasmonics 
or  metal  optics.  A  particularly  interesting  area  is  that  of 
plasmonics  in  the  terahertz  (THz)  regime,  where,  e.g.,  re¬ 
cent  research  on  the  excitation  of  surface  plasmon  modes 
in  tapered  metal  wire  geometries5-7  demonstrated  imaging 
of  electronic  characteristics  of  materials  on  micro-  and 
nanoscales.7,8 

Improvement  of  the  time-resolved  nanoscale  imag¬ 
ing  at  THz  frequencies  requires  an  efficient  local  emit¬ 
ter/detector  combined  with  broadband  time-domain  spec¬ 
troscopy  methods.9,10  Nano  structured  semiconductors,  in 
combination  with  recent  techniques  for  high-power,  ultra¬ 
fast  THz  excitation  sources,11-14  can  serve  as  intense  THz 
nanoscale  emitters.  Nanowires  have  been  the  subject  of 
extensive  ultrafast  spectroscopic  studies  at  visible  and  near- 
infrared  frequencies.15  Recent  experiments  characterized  the 
far-infrared  optical  response  of  GaAs  nanowires  using  ultra¬ 
fast  THz  spectroscopy.16  The  emission  of  THz  pulses  has 
been  reported  from  a  variety  of  nanostructures,17-23  but  the 
nature  of  the  ultrafast  excitation  mechanism  is  still  not  well 
understood.17-20  In  this  work,  we  show  that  InAs  nanowires 
can  serve  as  highly  efficient  radiators  of  THz  frequency 
electromagnetic  transients.  We  argue  that  nano  wires  overcome 
the  problem  of  internal  reflection  that  occurs  at  the  air- 
semiconductor  macroscale  interface,  leading  to  much  higher 
emission  efficiency  compared  to  planar  substrates.  Analysis 
of  the  dielectric  function,  including  transverse  confinement, 


suggests  that  the  emission  originates  from  a  low-energy  surface 
plasmon-polariton . 

II.  NANOWIRE  SAMPLES 

InAs  nanowires  are  grown  by  metal-organic  vapor  phase 
epitaxy  using  the  vapor-liquid-solid  (VLS)  technique.  The 
VLS  growth  mechanism  for  InAs  nanowires  has  been  thor¬ 
oughly  investigated.24,25  Our  version  of  this  process  uses  a 
Si-doped,  /2-type  GaAs  wafer  with  a  (111)  surface  orientation 
as  the  starting  substrate.  A  1-nm-thick  gold  deposition  is  made 
on  the  (lll)B  surface  by  electron  beam  evaporation.  Upon 
annealing  under  ASH3  at  600  °C,  the  gold  deposit  coalesces  to 
form  nanoparticles  that  represent  a  collector  for  InAs  growth 
precursors.  For  nanowire  growth,  the  temperature  is  reduced 
to  400  °C  and  trimethyl  indium  (TMIn)  is  introduced  into 
the  chamber.  The  nanowires  orient  vertically  as  shown  in 
the  scanning  electron  microscopy  (SEM)  image  in  Fig.  1(a). 
Their  length  is  in  the  range  10-20  pm,  the  diameter  tapers 
from  450  nm  at  the  base  to  60  nm  at  the  tip,  and  the  surface 
coverage  is  ~0.2  pm~2.  A  gold  nanoparticle  remains  attached 
to  the  tip  of  the  nanowire  after  growth.  In  addition  to  nanowire 
growth,  a  planar  wetting  layer  forms  due  to  competitive 
vapor-solid  heteroepitaxy  on  the  GaAs(lll)B  surface,  as 
shown  in  Fig.  1(b).  Wetting  layer  growth  occurs  at  a  rate 
smaller  than  the  vertical  nanowire  growth  rate,  leading  to  a 
thin  film.  We  estimate  a  maximum  wetting  layer  thickness 
of  200  nm  if  all  TMIn  is  consumed  in  this  pathway.  The 
presence  of  nanowires  implies  a  wetting  layer  thickness 
<  200  nm. 

III.  TERAHERTZ  SPECTROSCOPY:  RESULTS  AND 
DISCUSSION 

We  excite  the  InAs  nanowires  using  ultrashort  pulses  from 
a  mode-locked  Ti: Sapphire  laser  oscillator  (average  pulse 
energy  5  nJ,  duration  60  fs)  with  the  wavelength  centered  at 
820  nm.  The  //-polarized  pulses  are  focused  on  the  samples 
(~400  nJ/cm2)  at  an  incidence  angle  of  45°  with  respect 
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FIG.  1.  (Color  online)  (a)  Low-magnification  cross-sectional  SEM  image  of  the  InAs  nanowire  ensemble,  highlighting  the  nanowires’  high 
aspect  ratio  and  their  taper.  The  scale  bar  is  3  /xm  across,  (b)  Higher-magnification  SEM  image  showing  the  base  dimension  and  existence  of  a 
planar  wetting  layer  formed  by  vapor-solid  epitaxy.  The  scale  bar  is  1  /xm  across,  (c)  Sketch  (not  to  scale)  of  the  three  regions  of  the  sample 
that  are  potential  sources  of  THz  radiation  transients.  Our  experiments  show  that  the  signal  from  the  InAs  nanowires  dominates  the  signal  from 
the  InAs  wetting  layer  or  the  Si-doped  GaAs  substrate. 


to  the  substrate  surface.  We  characterize  the  emitted  electro¬ 
magnetic  transients  using  a  linear  autocorrelation,  where  two 
identical  pulses  are  delivered  to  the  sample  in  a  Michelson 
interferometer  geometry.26-28  The  THz  radiation  is  collected 
in  the  specular  direction  and  imaged  onto  a  liquid-He  cooled 
Si-bolometer  with  a  pair  of  off-axis  parabolic  mirrors.  A  Teflon 
filter  is  inserted  into  the  THz  beam  path  (unpurged  laboratory 
air)  to  remove  background  laser  light.  Linear  autocorrelations 
and  the  corresponding  power  spectra  (Fourier  transforms)  are 
presented  in  Fig.  2.  The  amplitude  of  the  signal  (at  zero  delay) 
obtained  with  the  nano  wires  (Fig.  2(a),  top)  is  only  about  two 
times  weaker  than  a  slightly  zz-doped  InAs  substrate  (Fig.  2(a), 
middle),  even  though  the  fill  factor  of  the  wires  is  only  of 
the  order  of  0.03.  This  shows  that  nanowires  are  unusually 
good  emitters  of  THz  radiation,  with  an  estimated  radiation 
enhancement  of  ~15  times  in  emitted  power,  compared  to  the 
InAs  substrate. 

To  verify  that  the  emitted  radiation  originates  from  the 
nano  wires,  we  test  a  control  sample  with  the  nanowires 
removed:  a  thin  (<200  nm)  InAs  wetting  layer  on  the  same 
GaAs  substrate  used  for  nanowire  growth.  To  remove  the 
nano  wires,  we  place  the  sample  in  an  ultrasonic  bath  and 


verify  their  absence  with  an  optical  microscope.  This  control 
sample  is  essential  to  isolate  the  contribution  of  the  InAs 
wires  from  the  InAs  wetting  layer,  which  builds  up  during  the 
growth  process.  The  control  sample  shows  a  very  weak  signal 
(Fig.  2(a),  bottom),  indicating  InAs  nano  wires  are  the  source  of 
strong  THz  radiation.  (The  emission  peak  occurring  ~3.5  THz 
demonstrates  the  broadband  response  of  the  detection  system.) 

To  understand  the  large  enhancement  of  the  THz  emission 
in  nanowires  compared  to  the  bulk  material,  we  must  first 
discuss  the  origin  of  the  THz  radiation.  This  has  two  com¬ 
ponents:  the  initial  excitation  mechanism  that  leads  to  carrier 
acceleration  and  the  subsequent  radiation.  In  bulk  semiconduc¬ 
tors,  ultrashort  optical  excitation  can  lead  to  dynamical  charge 
transport  of  photocarriers  that  drives  emission  of  electromag¬ 
netic  radiation.26,29  Excitation  of  the  semiconductor  above  the 
direct  energy  band  gap  generates  electron-hole  pairs  at  a  pene¬ 
tration  depth  determined  by  the  inverse  absorption  coefficient. 
In  the  case  of  InAs  nano  wires,  we  consider  four  physical 
mechanisms:  (1)  acceleration  in  a  surface  depletion  field; 
(2)  acceleration  due  to  fields  at  the  gold-nanoparticle/nanowire 
interface;  (3)  ambipolar  diffusion,  i.e.,  the  photo-Dember 
effect;29,30  and  (4)  nonlinear  mechanisms  such  as  optical 


delay  (ps)  frequency  (THz) 


FIG.  2.  (Color  online)  Autocorrelation  sig¬ 
nals  (a)  and  corresponding  Fourier  spectra 
(b)  for  an  InAs  nano  wire  sample  (top),  an  InAs 
substrate  (middle),  and  control  sample  (bottom), 
as  described  in  the  text.  The  top  two  signals  are 
plotted  on  the  same  relative  scale;  the  bottom 
signals  are  magnified  50  times.  Curves  are  offset 
vertically  for  clarity. 
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FIG.  3.  (Color  online)  THz  power  changes  linearly  as  excitation 
intensity  is  varied  (via  a  variable  neutral  density  filter)  over  a  range  of 
150.  The  inset  shows  the  normalized  spectra  of  THz  transients  from 
InAs  nano  wires  as  the  intensity  varies  by  a  factor  of  25.  Absence  of 
a  spectral  shift  is  consistent  with  radiation  emanating  from  coherent 
plasma. 

rectification  (OR)  or  shift  currents.31  The  first  mechanism  is 
known  to  be  important  in  ft-doped  GaAs  and  InP,  where  the 
surface  depletion  field  is  large  (>10  kV/cm). 32,33  Narrow-gap 
semiconductors  such  as  Te,  InSb,  and  InAs  have  negligible 
surface  fields,  making  this  mechanism  unimportant  here.30 
The  second  mechanism  is  also  insignificant  because  the  Fermi 
level  is  pinned  at  the  conduction  band  minimum  in  InAs, 
resulting  in  a  negligible  Schottky  barrier  and  band  bending.34 
If  a  Schottky  barrier  did  exist,  the  associated  depletion  region 
would  only  account  for  a  tiny  volume  fraction  of  each 
nano  wire. 

To  differentiate  between  the  photo-Dember  and  the  nonlin¬ 
ear  processes,  we  note  that  the  power  of  THz  pulses  originating 
from  OR  or  shift  currents  would  increase  quadratically  with 
laser  excitation  at  small  conversion  efficiencies  and  saturate  at 
high  intensities.  The  power  dependence  of  the  photo-Dember 
effect  can  be  complicated  because  transport  depends  critically 
on  carrier  density  through  scattering  and  screening.  For 
bulk  systems  and  under  various  excitation  conditions,  the 
dependence  of  the  radiated  THz  power  on  the  excitation  levels 
has  been  observed  to  be  either  linear29,35  or  quadratic36,37 
at  low  power,  while  exhibiting  saturation  behavior  at  higher 
powers.35-37  Figure  3  shows  that  the  emitted  power  increases 
linearly  with  excitation  intensity  as  the  pumping  level  varies 
by  more  than  two  orders  of  magnitude.  This  rules  out  any 
meaningful  contribution  from  OR  or  shift  currents,  leaving 
only  the  photo-Dember  mechanism  as  a  plausible  source  for 
the  initial  carrier  acceleration.  We  also  point  out  that  the 
observed  emission  was  largely  independent  of  the  polarization 
state  of  the  excitation,  in  accord  with  the  photo-Dember 
mechanism  and  in  contrast  to  that  expected  from  the  nonlinear 
contributions. 

While  the  photo-Dember  effect  can  directly  lead  to  radia¬ 
tion,  it  is  also  well  known  that  charge  carriers  will  move  to 
screen  the  photocarrier  distribution.  If  the  perturbation  occurs 
on  a  time  scale  shorter  than  the  inverse  cold  plasma  frequency 
and  the  rate  of  momentum  relaxing  collisions  is  sufficiently 
small,  coherent  collective  motion  of  cold  charge  carriers 
(plasmons)  will  take  place,  leading  to  a  macroscopic  dipole 
radiating  close  to  the  plasma  frequency.26,38  Experiments  with 
bulk  InAs  and  GaAs  planar  layers  identified  cold  plasmon, 


phonon,  and  mixed-mode  components  in  the  THz  radiation 
spectra.26,28,39  These  longitudinal  modes  are  poles  of  the 
polarization  and  must  be  included  in  the  analysis  of  the 
free- space  radiation  spectra,  regardless  of  the  exact  nature  of 
the  excitation  mechanism.4 

We  assume  that  the  radiation  originates  from  plasmons  and 
use  a  classical  Drude  response38  to  analyze  the  spectra  in  Fig.  2. 
The  simplest  approach  is  to  first  consider  the  case  of  bulk 
plasmons.  The  plasmon  is  excited  at  frequency  co  =^/ o)2p-y 2/4, 
where  y  is  the  damping  rate  (i.e.,  momentum  relaxation  time) 
attributed  to  electrons  scattering  with  longitudinal-optical 
phonons  and  cop  is  the  Langmuir  plasma  frequency: 


n.e* 


m*£o£oo 


(1) 


Here,  ne  is  the  electron  density,  So  is  the  permittivity  of 
free  space,  m*  =  0.023me  is  the  electron  effective  mass  in 
the  conduction  band,  and  =  11.8  for  InAs.  Assuming 
an  electron-donor  ion  plasma  (i.e.,  neglecting  holes),  we 
estimate  an  electron  density  for  the  wire  ensemble  to  be 
~5  x  10 15  cm-3  based  on  the  frequency  of  the  peak  maximum 
and  its  broadening  [Fig.  2(b)]. 


IV.  ELECTRONIC  TRANSPORT  MEASUREMENTS 


To  test  whether  the  bulk  plasmon  theory  explains  our 
results,  we  independently  measure  the  carrier  concentration 
by  performing  electronic  transport  measurements.  This  is 
accomplished  with  two  independent  methods:  (1)  probing 
individual  nanowires  directly  on  the  growth  substrate  inside  of 
a  SEM41,42  and  (2)  fabricating  individual  nanowire  field  effect 
transistor  (FET)  devices.  Individual  nanowires  are  contacted 
with  a  tungsten  nanoprobe  retrofitted  inside  of  a  SEM.  A  large 
area,  backside  ohmic  contact  completes  the  electric  circuit 
[Fig.  4(a)].  The  high  aspect  ratio  geometry  of  the  nanowires 
leads  to  reduced  carrier  screening  and  the  onset  of  a  space 
charge-limited  current  (SCLC)  at  lower  voltage  than  would  be 
expected  for  a  bulk,  low-aspect  ratio  specimen  with  a  similar 
free  carrier  concentration  43  The  symmetric  nature  of  the  l-V 
curve  in  Fig.  4(b)  afforded  by  the  two  ohmic41,42  contacts 
despite  their  geometrically  asymmetric  nature,  is  consistent 
with  SCLC  (a  bulk- limited  regime).  This  is  in  marked  contrast 
to  rectifying  I-V  characteristics  observed  for  Au  catalyst/Ge 
nanowire  contacts,  where  a  high  Schottky  barrier  is  formed  44 
Transport  is  diffusive  in  SCLC,  so  for  a  specimen  with  a 
nanowire  geometry,  the  mobility  can  be  calculated  directly 
using  the  relation43 


fi  = 


LI  SCLC 

jtsV2  ’ 


(2) 


where  L  is  the  nanowire  length,  Isclc  is  the  current,  and 
V  is  the  applied  bias.  Once  the  mobility  is  known,  the 
carrier  concentration  can  be  extracted  from  the  conductivity 
(measured  at  low  bias). 

Alternatively,  ne  can  be  determined  using  the  crossover 
voltage  Vc  between  the  ohmic  and  SCLC  regimes.  Accounting 
for  SCLC  geometric  scaling  in  nanowires43,45  we  have 

ne  =  sVc/eR2,  (3) 


where  R  is  the  nanowire  radius. 
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FIG.  4.  (Color  online)  (a)  SEM  im¬ 
age  depicting  the  probing  of  an  individual 
nanowire  on  the  growth  substrate.  The 
scale  bar  corresponds  to  1.5  /xm.  (b)  I-V 
curve  for  the  InAs  nanowire  measured  in 
(a).  In  the  inset,  the  same  data  plotted  in 
the  formalism  of  the  SCLC  enable  us  to 
extract  the  ohmic  (intercept)  and  SCLC 
(slope)  contributions  to  the  total  current 
and  hence  the  crossover  voltage. 


In  the  second  experiment,  source  and  drain  Ti/Au  electrodes 
are  fabricated  over  randomly  dispersed  nanowires  on  an 
n- Si/Si02  substrate.  An  example  of  one  of  the  nanowires, 
along  with  its  electric  transport  data,  is  shown  in  Fig.  5. 
We  calculate  the  coupling  capacitance  C  between  the  InAs 
nano  wire  and  the  back  gate  through  the  100-nm-thick  SiC>2 
gate  dielectric  using  the  cylinder-on-plate  model.  Mobility  and 
carrier  concentrations  are  then  extracted  from  the  transcon¬ 
ductance  using  the  relations  /x  =  p— and  ne  =  CV±/e 
where  Vsd>  Vg,  and  Vth  are  the  source  drain,  gate,  and  threshold 
voltage,  respectively.  Both  measurement  techniques  show  that 
as  the  nanowire  diameter  increases  from  ~40  to  ~250  nm,  the 
carrier  concentration  decreases  from  ~1018  to  ~1016/cm3  and 
the  mobility  increases  from  ~100  to  3000  cm2/V-s  41  These 
trends  are  in  reasonable  agreement  with  other  measurements 
on  VLS -grown  InAs  nanowires  using  a  FET  geometry.46 
However,  nanowire  mobility  measurements  provide  lower- 
limit  estimates  of  the  intrinsic  mobility,  because  ohmic  losses 
at  the  contacts,  as  well  as  surface  defects  introduced  during 
processing  of  the  nanowires,  tend  to  lower  the  measured 
values. 


V.  PLASMON  DISPERSION  IN  NANOWIRES 

Our  electronic  measurement  of  the  carrier  concentration 
differs  by  nearly  two  orders  of  magnitude  from  that  determined 
by  the  bulk  model  Eq.  (1),  suggesting  that  the  bulk  model 
is  not  applicable  to  the  nanowire  geometry.  To  resolve  this 
discrepancy,  we  consider  a  more  realistic  description  of  the 
longitudinal  polarization  modes  for  the  cylindrical  nanowire 
geometry  by  applying  a  formalism  developed  by  Pitarke  et  al.41 
This  theory  provides  the  dielectric  function  of  a  conducting 
cylinder,  yielding  eigenmodes  of  the  geometric  structure  as 
a  function  of  parameter  q^R  and  an  integer  mode  number 


m ,  arising  from  the  cylindrical  symmetry.  The  dimensionless 
parameter  q^R  is  the  product  of  the  wavevector  component 
along  the  nano  wire  axis  qi  and  the  nanowire  radius  R.  We 
modify  the  model  to  consider  discrete  axial  modes  with  qi  = 
in /L,  where  L  is  the  length  of  the  nanowire  and  i  >  0  is 
an  integer 48,49  Figure  6(a)  shows  the  dispersion  relation  of 
the  modes  for  the  case  of  the  assumed  local  response.47  In 
addition  to  the  dispersionless  bulk  mode  at  cop,  there  exist 
surface  plasmon  modes  indexed  by  the  integer  m  and  given 
by47’50 


co(m,qiR)  =  cop 


&OG 


(4) 


where  am(x)  =  I'm(x)Km(x)/ Krm(x)Im(x),  Im  and  Km  are 
modified  Bessel  functions,  and  the  prime  indicates  the  deriva¬ 
tive  of  the  corresponding  functions.  In  the  limit  of  q^R^> 
1  (am  — >  —  1),  the  surface  plasmon  modes  approach  the 
planar  surface  plasmon  frequency  &>sp  =  ^  V£oo/(£oo  +  £o) 
[Fig.  6(a)],  which  reduces  to  the  familiar  <x>sp  =  cop/V2  for 
£0  =  =  1 .  The  m  =  0  modes  exhibit  acousticlike  dispersion 

for  small  q^R,  while  the  higher  modes  (m  >  0)  are  essentially 
dispersionless  around  the  cosp  frequency.  Thus,  for  the  assumed 
purely  longitudinal  excitation,  the  only  excitable  radially 
symmetric  m  =  0  mode  is  considered  in  the  remainder  of 
our  analysis. 

To  analyze  the  radiated  spectra  of  the  plasmon  modes, 
we  consider  the  induced  polarization  by  an  externally  ap¬ 
plied  longitudinal  electric  field  Eext  (such  as  photo-Dember) 
given  by49 


P(q,co)  =  £oo(l - - — - 

V  e(q,<o) 


Eext(q,co), 


(5) 
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FIG.  5.  (Color  online)  (a)  Source-drain  cur¬ 
rent  vs  source-drain  voltage  sweeps  of  a  rep¬ 
resentative  FET  at  different  fixed  gate  biases 
showing  the  n- type  depletion  mode  behavior, 
(b)  Transfer  characteristic  of  the  FET  shown  in 
(a),  which  enables  determination  of  the  mobility 
and  carrier  concentration  by  way  of  the  transcon¬ 
ductance  and  threshold  voltage,  respectively.  The 
extracted  device  parameters  are  ne  =  4.8  x 
1017  cm-3  and  /x  =  101  cm2/V-s.  Nanowire 
dimensions  are  R  =  34  nm  and  L  =  3  /.xm. 
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q*R 


frequency  (THz) 


FIG.  6.  (Color  online)  (a)  Dispersion  of  the 
longitudinal  bulk  and  surface  plasmon  modes, 
(b)  Normalized  measured  emission  spectrum  of 
the  InAs  nanowire  ensemble,  (c)  Theoretical 
spectra  of  the  lowest-order  axial  surface  acoustic 
plasmon  modes  of  the  nanowire,  with  ne=  4.8  x 
1017  cm-3,  R  =  34  nm,  L  =  3  /xm,  and  plasmon 
broadening  of  0.05 cop. 


where  £(q,co)  is  the  wavevector  and  frequency-dependent 
dielectric  function  of  the  wires.47  The  radiated  far  field 
is  £ra d  oc  d2P/dt2,  so  ignoring  propagation  effects  and  as¬ 
suming  broadband  excitation,  the  radiated  power  spectrum 
Srad(g,&0  °C  |1  —  s~x{q,co)\2coA .  Fig.  6(c)  shows  two  such 
spectra  for  axially  symmetric  i  —  1,3  modes,  taking  R/L  ~ 
0.01  of  the  nano  wire  shown  in  Fig.  5.  In  accordance  with 
antenna  theory,  only  the  axially  symmetric  modes  (I  =  odd) 
interfere  constructively  in  the  far  field.48,49  The  frequency 
axis  in  Fig.  6(c)  is  scaled  by  the  bulk  plasmon  frequency 
of  ~12  THz  as  determined  from  Eq.  (1)  for  ne  =4.8  x 
1017  cm-3.  The  predicted  acoustic-plasmon  frequency  of  the 
(m  =  0,  i  =  1)  mode  in  this  nanowire  geometry  is  within  a 
factor  of  two  of  the  experimentally  observed  value  [Fig.  6(b)]. 
The  next-lowest  axial  mode  shows  even  better  overlap  with 
the  measured  spectrum.  This  is  in  stark  contrast  to  the 
factor-of-100  discrepancy  obtained  from  the  bulk-plasmon 
model  analysis  Eq.  (1)  and  suggests  that  the  observed  THz 
emission  originates  from  a  low-energy  acoustic-plasmon  mode 
of  the  nanowires.  This  interpretation  is  also  consistent  with 
independent  measurements  of  electric  transport  and  ultrafast 
THz  spectroscopy. 

We  point  out  that  the  experimentally  observed  THz  spectral 
shape  is  inhomogeneously  broadened  due  to  the  ensemble 
of  nanowires  of  varying  geometry  and  carrier  concentra¬ 
tion.  Variations  in  the  aspect  ratio  ( R/L )  would  result  in 
a  proportional  linear  shift  of  the  frequency  of  the  acoustic 
modes  [Fig.  6(a)],  while  the  carrier  density  variations  would 
correspond  to  quadratic  frequency  shifts.  Averaging  over  the 
nanowire  ensemble  smears  out  the  spectrum,  and  prevents 
observation  of  distinct  features  associated  with  an  individual 
nanowire,  as  shown  in  Fig.  6(c).  We  expect  the  tapering  of  the 
nano  wire  to  become  important  when  q  A  R  1 .  In  our  system, 
the  lowest-order  mode  has  qAR  ~  0.03  1,  suggesting  that 

the  approximation  of  a  constant  radius  is  appropriate. 

Our  experimental  results  and  modeling  suggest  that  the 
presence  of  a  low-energy  acoustic  surface  plasmon  mode  is 
responsible  for  THz  emission  from  InAs  nanowires.  Surface 
plasmon  cavity  modes  in  metallic  nanowires  have  been  re¬ 
cently  inferred  by  the  near-field  microscopy51  and  optical  loss 
spectroscopy.48  Our  data  indicate  that  we  observe  radiation 


directly  from  the  acoustic  surface  plasmon  modes  in  the 
Hertzian  dipole  geometry.  These  results  are  corroborated  by  a 
recent  observation  of  THz  emission  from  Si  nanowires,  where 
the  role  of  surface  modes  also  was  discussed.20  In  addition, 
we  were  not  able  to  detect  any  THz  emission  from  much 
shorter  (L  ~  300  nm)  nanowires,18  where  an  appreciable 
momentum  mismatch  for  radiation  outcoupling  cannot  be 
compensated  by  the  lowest  axial  modes,  even  in  the  Hertzian 
dipole  approximation.  Efficient  Landau  damping  of  plasmons 
is  also  likely  to  exist  in  short  wires,  further  reducing  THz 
emission.52 


VI.  RADIATION  OUTCOUPLING 

Finally,  we  address  the  observed  large  enhancement  of 
the  emitted  radiation  from  nanowires  in  comparison  to  a 
planar  substrate.  For  an  InAs  substrate,  coherent  plasmons 
are  excited  below  the  surface,  leading  to  charge  oscillations 
in  the  direction  of  the  surface  normal  and  a  concomitant 
emission  of  radiation  into  free  space.  The  sin2(0)  dipole 
radiation  pattern  is  sketched  in  Fig.  7(a),  where  0  is  the  internal 
angle  with  respect  to  the  surface  normal.  The  radiation  pattern 
is  primarily  oriented  parallel  to  the  semiconductor  surface. 
When  propagating  from  a  high  refractive  index  semiconductor 
(n  ~  4)  into  free  space,  Snell’s  law  at  the  interface  must  be 
satisfied  or  total  internal  reflection  will  occur.  Outcoupling 


ray  1 


InAs  substrate 


f  nanowire 


InAs  wetting  layer _ 

GaAs'(1 1 1’Js'substrat^™ 


FIG.  7.  (Color  online)  (a)  Radiation  pattern  from  a  plasma  dipole 
(oscillation  direction  indicated  by  bold  double  arrow)  beneath  the 
surface  of  a  planar  semiconductor.  The  shaded  region  within  the  lobes 
depicts  the  radiation  escape  cone;  only  light  within  this  cone  (ray  1) 
can  escape  into  free  space.  The  remaining  radiation  (ray  2)  is  trapped 
due  to  total  internal  reflection,  (b)  Schematic  of  coherent  plasmon 
motion  and  light  propagation  with  a  freestanding  semiconductor 
nanowire. 
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from  the  dipole  takes  place  in  a  narrow  escape  cone  with  an 
apex  angle  of  0C  =  arcsin(l/n)  ~  15°,  sketched  as  the  shaded 
triangle  in  Fig.  7(a).  Light  rays  striking  the  surface  at  angles 
6  greater  than  the  apex  angle  of  the  escape  cone  suffer  total 
internal  reflection  and  remain  trapped  in  the  semiconductor. 
Only  a  small  portion  of  the  dipole  radiation  pattern  (<1%) 
has  the  proper  wavevector  to  satisfy  coupling  to  free  space.30 
The  refractive  index  is  a  complex  quantity  and  can  be  highly 
dispersive  in  the  vicinity  of  the  plasma  and  optical  phonons. 
Plasma  oscillations  are  strongly  damped  in  In  As,  however,  so 
the  complex  index  is  relatively  constant  at  and  above  1  THz. 

The  radiation  trapping  problem  is  well  understood  and  has 
led  to  the  development  of  different  schemes  for  improving 
the  extraction  efficiency.  For  example,  the  direction  of  charge 
carrier  motion  has  been  altered  by  the  Lorentz  force,53  where 
a  large  constant  magnetic  field  (~1  T)  rotates  the  dipole 
to  preferentially  orient  one  lobe  of  the  radiation  pattern 
in  the  escape  cone.30,53  Another  approach  is  to  implement 
structured  surfaces,  such  as  low  absorption,  index-matched 
prism  stripes54  or  metal  edges.55  In  contrast,  freestanding, 
vertically  oriented  nanowires  possess  an  inherent  geometric 
efficiency  compared  to  substrates  because  their  primary  sur¬ 
face  is  parallel  to  the  direction  of  collective  charge  motion.  The 
effect  is  sketched  in  Fig.  7(b).  We  estimate  the  enhancement  of 
nanowire  emission  over  the  substrate  by  calculating  the  ratio  of 
emitted  power  into  the  full  solid  angle  from  nanowires  (Pnw) 
to  that  emitted  (and  coupled  to  free  space)  from  the  substrate 
(PBuik): 

^nw  _  f  dQ  sin 2(0) 

P^k~  J dQsmH0)(l  -  Rll(0)y  U 

where  the  assumed  source  is  a  Hertzian  dipole  emitter 
(a  sin2(0))  and  dQ  =  sin  BdBdtp  is  an  infinitesimal  solid 
angle  in  spherical  geometry.  We  account  for  the  standard 
Fresnel  losses  (1  —  R\\{0))56  when  considering  emission  from 
photo-Dember  excitation  (i.e.,  p-polarized  by  the  nature  of 
dipole  orientation)  in  the  substrate  (denominator  of  Eq.  (6)). 
Because  the  nano  wire  volume  is  much  smaller  than  the  cube 
of  the  characteristic  THz  wavelength,  its  effective  refractive 
index  is  near  unity;  hence,  no  Fresnel  losses  are  considered 
(numerator).  The  integration  over  0  in  the  case  of  a  substrate 
is  performed  from  0  to  0C  to  account  only  for  the  escaped 
emission.  Using  Eq.  (6)  for  index  n  ~  4,  the  estimated 
enhancement  factor  is  ~20  and  is  practically  unchanged,  even 
when  the  experimental  collection  geometry  (45°  incidence 
with  10°  collection  half-angle)  is  taken  into  the  account.  The 


observed  emission  enhancement  per  nanowire  of  a  factor  of 
15  compares  favorably  with  this  estimate. 

An  enhancement  of  THz  emission  has  been  observed  in 
several  other  material  systems.  A  recent  experiment  with 
opaque/blackened  silicon  showed  that  surface  roughness  in 
the  form  of  needles  (length  of  ~2  /zm,  diameter  of  ~300  nm) 
enhanced  the  THz  output  compared  to  an  untreated  surface.19 
This  was  ascribed  to  increased  absorption  of  pump  laser 
photons  via  multipath  scattering.  The  emission  mechanism 
was  attributed  to  an  acceleration  of  photocarriers  in  the  surface 
depletion  field  region.19  Radiation  power  from  blackened 
silicon,  however,  was  ~50  times  weaker  than  a  polished 
InAs  substrate.  THz  enhancement  was  also  observed  with  InN 
nanopillars  for  a  surface  number  density  that  is  ~250  times 
larger  than  in  the  present  work.17  Our  experiments  show 
that  InAs  nanowires  generate  signal  levels  comparable  to 
substrates,  not  even  accounting  for  the  spatial  fill  factor. 
This  performance  makes  InAs  nanowires  highly  attractive  for 
nanoscale  THz  emitter  and  detector  applications. 

In  summary,  we  have  demonstrated  strong  emission  of  THz 
radiation  from  a  cold  plasma  in  InAs  nanowires.  There  is 
a  15 -fold  enhancement  in  power  of  outcoupling  efficiency 
compared  to  a  bulk  InAs  planar  substrate,  which  we  attribute 
to  a  more  favorable  geometry.  A  geometry-dependent  dielec¬ 
tric  function  that  includes  a  low-frequency  acoustic  surface 
plasmon  mode  provides  a  good  description  of  the  observed 
THz  radiation  spectra  and  order  of  magnitude  agreement  with 
separate  electronic  measurements.  A  systematic  study  of  the 
THz  emission  from  single  nanowires  as  a  function  of  geometry 
must  be  performed  to  better  understand  plasmon  dispersion  in 
these  structures. 
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